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Abstract. We investigate the scaling relations of bulge and disk structural parameters for a sample of 108 disk galaxies. 
Structural parameters of individual galaxies are obtained from two-dimensional bulge/disk decomposition of their H-band 
surface brightness distributions** . Bulges are modelled with a generalized exponential (Sersic) with variable integer shape index 
n. We find that bulge elfective scalelengths rf and luminosity increase with increasing n, but disk properties are independent 
of bulge shape. As Hubble type T increases, bulges become less luminous and their mean effective surface brightness < /i^ > 
gets fainter; disk < yuf > shows a similar, but much weaker, trend. When bulge parameters (< yuf >, , M®) are compared with 
disk ones (< /jf >, rf , M"), they are tightly correlated for n = 1 bulges. The correlations gradually worsen with increasing n 
such that n = 4 bulges appear virtually independent of their disks. The Kormendy relation, < yuf > vs. rf , is shown to depend 
on bulge shape n\ the two parameters are tightly correlated in n = 4 bulges (r = 0.8), and increasingly less so as n decreases; 
disk < fi^ > and rf are well correlated (r = 0.7). Bulge-to-disk size ratios /rf are independent of Hubble type, but smaller 
for exponential bulges; the mean rf /rf for n = I bulges is 4 times smaller than that for n = 4, with a spread which is 9 times 
smaller. Strongly barred SB galaxies with exponential bulges are more luminous than their unbarred counterparts. Exponential 
bulges appear to be closely related to their underlying disks, while bulges with higher n values are less so; n = 4 bulges and 
their disks apparently have no relation. We interpret our results as being most consistent with a secular evolutionary scenario, 
in which dissipative processes in the disk are responsible for building up the bulges in most spirals. 

Key words, galaxies: evolution - galaxies: formation - galaxies: fundamental parameters - galaxies: spiral - galaxies: structure 
- infrared:galaxies 

1 . Introduction with measured rotation curves, and //) postponed to future work 

the discussion of the impact of scaling relations of bulges and 

This is the fifth paper in a series based on near-infrared imaging ■■ , »u ■ * j w t »u * 

" ° disks on their formation and evolution. In the present paper, we 

of a sample of disk galaxies m the Pisces-Perseus supercluster. n jj .u - • 

^ ° ,_L . will address this issue. 

Akeady in Paper III of the series (Moriondo et al. 1999b i, we 

explored the relations between structural parameters of bulges The traditional picture is that bulges formed well before 
and disks, but in the context of distance measurements and the the disk of a galaxy, in a "monolithic collapse", either as the 
consequent definition of their Fundamental Planes. For this rea- first rapid stage of the collapse of a galaxy-size density per- 
son, that paper: /) made only use of the subset of galaxies (40) turbation, as in the Milky Way collapse model of Eggen et 



al. (I1962> . or through the infall of enriched gas from the star- 

Send offprint requests to: L. K. Hunt ^^^.j^jj^g j^^j^ ^^^^^ .g^.jy„ ^^^j^^ forj^ation is 

* Based on observations at the TIRGO, NOT, and VATT telescopes. , j- ^ j -^u- r i i_- i u j i 

,^ . ,, ' „ . ' . . ^. also predicted within the framework of hierarchical cold-dark 

TIRGO (Gornergrat, CH) is operated by IRA-CNR, Arcetri, Firenze. ^ ,^ ■ r ^ ■ , , • 

NOT (La Palma, Canary Islands) is operated by NOTSA, the Nordic "^^"^i" ^^DM) scenarios of structure formation and evolution 

Observatory Scientific Association. VATT (Mt. Graham, AZ) is oper- (KauflFmann et al. EMI Baugh et al. 1996 ). In this case, the 

ated by VORG, the Vatican Observatory Research Group. bulge originates from the merging of stellar disks within the 

** Table 1 is only available in electronic form at the CDS merging halos, while the disk that is observed today is assem- 

via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via bled through gas accretion at later epochs. A more recent vari- 

thttp://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-)-A/, ation of this scenario is the growth of galactic bulges through 
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merging of dense satellites (Aguerri et al. '200111. Evidence in 
favor of "early" bulge formation rises from studies of the stellar 
populations in the Milky- Way bulge (Ortolani et al. 1995 1 and 
in bulges of nearby galaxies (Jablonka et al. 119961. In addition, 
bulge velocity fields (Kormendy & lllingworth 1982) and their 
high central surface brightness seem to be imprints of a forma- 
tion process that involved dissipational collapse (e.g., Carlberg 
I1986> at epochs earlier than the formation of the disks and, per- 
haps, of the stellar halo. 

More recently, it has been suggested that the bulges in 
(some) disk galaxies can be formed after the disk, through sec- 
ular evolutionary phenomena. Numerical simulations (Combes 
& Sanders .1981. Combes et al. .l990.. Norman et al. 1996 1 have 
shown that a spheroidal stellar component can originate from 
dissipationless processes such as the thickening or destruction 
of a bar Moreover, when gas is included in the simulations, dis- 
sipation accelerates the evolution, giving rise to galaxies which 
are more centrally concentrated (Junqueira & Combes 1996 1. 
Indeed, some of the observed bulge characteristics such as tri- 
axial or peanut shape and cylindrical rotation (e.g.. Shaw l 19^ 
Liitticke et al. 2000 1 are predicted by such models. However, 
gas is not necessary for dissipation, as pointed out by Zhang 
in a series of papers (Zhang [T996I [T9981 [T999t . Secular evo- 
lution may also operate through outward transport of angular 
momentum through a collective dissipative process linked to 
spiral density waves. Observational evidence in favor of the 
secular scenario include the similarity of bulge and disk color 
indices (e.g., Balcells & Peletier 1994,1, the mass-to-luminosity 
ratios of the components (Widrow et al. I2003> . the disk-like 
kinematics observed in many bulges (Kormendy 1993 1, and 
the correlation of bulge and disk scale lengths (de Jong 1996 1. 
Indeed, when a simple exponential is used to fit the bulge sur- 
face brightness profile, a strong correlation is found between 
bulge and disk scalelengths (Courteau et al. 1996 1, seemingly 
supporting the secular formation scenario for bulges. The cor- 
relation was later shown to depend critically on the use of a 
simple exponential as a bulge fitting function, since when gen- 
eraUzed exponentials were used instead, the mean scalelength 
ratio depends on the bulge shape index (Moriondo et al. 1998a| 
- hereafter MGH; Graham 20Q1 1. Thus it appears, albeit for 
small samples, that the ratio of bulge and disk scalelengths 
may not be constant over a range of Hubble types as previ- 
ously claimed. This is of importance not only for the chronol- 
ogy of bulge and disk formation, but also for the meaning of 
the Hubble sequence itself in terms of photometric observables 
(scale-free Hubble sequence). 

To further investigate bulge and disk structural parameters, 
and to better constrain their connection with formation scenar- 
ios, we have performed two-dimensional (2D) bulge/disk de- 
compositions on a large sample of disk galaxies in the Pisces 
Perseus supercluster (see Moriondo et al. I1999al hereafter 
Paper I). These decompositions are performed on image data 
in the H band because the near-infrared (NIR) wavelengths are 
the most direct tracer of the stellar mass (e.g., Rix & Rieke 



I1993i Gavazzi I1993I I. Moreover, NIR emission is much less 
affected by recent episodes of star formation which typically 
dominate disk morphology at optical wavelengths. Finally, the 
effects of dust extinction are mitigated at 1 .65 //m (H) since Ah 



is roughly a factor of 8 lower than Ay- We describe the sample 
and the imaging data in Sect.|2] together with the bulge/disk de- 
compositions. Section 13 presents our results, and we interpret 
them in Sect.|?]in the context of various evolutionary scenarios. 

2. The sample, the images, and the structural 
decomposition 

The galaxies under investigation belong to a sample selected 
from the catalogued disk galaxies in the Pisces-Perseus (PP) 
supercluster area, 22*" < RA < 4^ and 0° < 6 < 45°, excluding 
obvious foreground members with heliocentric velocities less 
than 3000 km s"'. The original sample was restricted to about 
950 galaxies with optical major axis 0'.5 < a < 4'.0 and avail- 
able good-quality 21 -cm spectra (Giovanelli & Havnes [19891 
and references therein). From this set two distinct subsamples 
were extracted. Sample A comprises about 150 galaxies, ran- 
domly selected to cover the entire area and inclination range, 
and with types Sb or later. Sample B comprises the 50 clus- 
ter galaxies within an angular distance of 5.5° from the Pisces 
cluster center (l** 20™, 33°), with optical size a > 0'.8 and in- 
clination 30° < i < 75°, spanning the range SO to Sd in Hubble 

type- 
Out of these 200 objects, images were obtained for 174 
galaxies (see Paper I) in the H band (35 in JHK) with the 
ARNICA camera (Lisi et al. 1993 1996 Hunt et al. [T996t 
mounted on various telescopes. We refer the reader to Paper 
I for details concerning observations, data reduction, and pho- 
tometric calibration. Analysis of the JHK color images is pre- 



sented in Moriondo et al. (120011 Paper IV). Unless otherwise 
specified, all the photometric data derived and here reported 
refer to the set of images in the H band. 

2.1. Two-dimensional bulge/disk decomposition 

As in MGH, we fitted the images with a two-dimensional (2D) 
parametric brightness distribution model, consisting of a gen- 
eralized exponential (Sersic 1968) bulge of shape index n : 



-a„ 



(x,y) - /<,exp 
and an exponential thin disk: 



(1-6*) 



B\2 
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, (1) 



7^(^,3') = /rf(0)exp 



1 

rd 



(2) 



/p and re are effective (half-light) values, is the apparent 
bulge ellipticity, a„ is a constant relating the effective bright- 
ness and radius to the exponential values (see MGH). x and y 
are in arbitrary units, with origin at the galaxy's center, and x 
along the major axis. The bulge is assumed to be an oblate rota- 
tional ellipsoid (Kent et al. ll991> . coaxial with the disk, and its 
apparent eccentricity is related to the intrinsic eccentricity 
e"^ by: 



= e'^sin; . 



(3) 
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Fig. 1. Upper panels: surface brightness profile along the major 
axis (from elliptical isophotal fitting) in H mag arcsec"^. Lower 
panels: radial run of the reduced;^'^ of the best-fitting 2D model 
for different values of the n. The first panel shows a galaxy with 
a best-fit n of 4, while a best-fit « = 1 is shown in the second 
panel. Abscissas are labelled in pixel, 1 pixel = 0.975 arcsec. 



The structural parameters have been determined by fitting 
the 2D image to the model, convolved with a circular Gaussian 
seeing disk of appropriate FWHM, using a minimization. 
The parameters left free to vary in the fit include: the two sur- 
face brightnesses Ig and ld{0) (or jje and yuo when given in mag- 
nitudes), the two scale lengths re and rd, the bulge ellipticity 



e , and the system inclination /. The fitting algorithm and tests 
of its reliability are described more fully in MGH. 

The use of a 2D fitting technique is mandatory in the anal- 
ysis of objects of middle to high inclination, such as samples - 
including the present one (68 galaxies with / > 40°) - selected 
for Tully-Fisher studies. A feature of our algorithm is its abil- 
ity to suppress non-axisymmetric structure. This is achieved 
through a "folding" procedure in which we exploit the model 
symmetry about the major and minor axes. The galaxy image 
is first "folded" about the major axis (measured a priori), in 
order to obtain the "long" semi-ellipse; it is then folded again 
about the minor axis in order to obtain only one quadrant. There 
are consequently four points used to infer the value of a sin- 
gle pixel in the folded quadrant, and from these, a mean and 
a standard deviation are calculated; only this single quadrant 
is input to the fitting routine. Typically pixels are weighted by 
the photon statistics, calculated by modelling the detector noise 
characteristics and assuming background limited performance 
(verified to be true in all cases). The weighting scheme how- 
ever gives greater weight to those regions with low standard 
deviation since it chooses the maximum of the Poisson statis- 
tics or the standard deviation resulting from the folding pro- 
cedure. The regions with significant asymmetry show up with 
large standard deviation in the folding process, and they are 
given accordingly less weight in the decomposition. 

Because of the known mutual functional dependency of r^, 
jde, and n, we chose to use a "grid method" for the determi- 
nation of n rather than include it as a free parameter (see also 
MacArthur et al. l2003l hereafter MACH03). Because of the 
large sample and consequent numerical calculation time for the 
2D fits, we restricted n to integer values of 1, 2, 3 or 4, and ran 
fits for each of these separately. The best-fitting n value was 
then derived by minimizing the reduced over the inner re- 
gion. Two examples of the radial run of xl for different n's are 
shown in Fig. [0 The surface brightness profiles in the upper 
panels are derived from elliptical isophotal fitting and are only 
used for display purposes; the lower panels show the behavior 
of the reduced resulting from the 2D best fit at different n, 
averaged within the same elliptical annuli of the upper panel. 
As can be seen, the bulge fit influences primarily the inner part 
of the galaxy. It is consequently necessary to consider the value 
of xl as a function of n especially at small radii, because the 
large area of the disk tends to make the influence of the bulge 
shape in the global x^ very small and easily masked by other 
features. The median (and mean) FWHM seeing of the images 
in our database is \'!9. Only 4 of 108 galaxies have rf smaller 
than this, but it is difficult to distinguish star clusters or nuclear 
point sources from high-n bulges. Point sources have not been 
included in our fits. 

The discretized range of integer n values 1, 2, 3, 4 sam- 
ples rather well the continuous range of n values found by 
previous authors. MACH03, on the basis of one-dimensional 
(ID) profile fits, conclude that the surface brightness profiles 
of Freeman Type-I galaxies (Freeman 1970 1 are adequately 
described by a double-exponential distribution with bulges of 
0.1 < « < 2; however only 8 of the 121 galaxies studied by 
them have n < 0.5. On the other hand, with 2D fits, MollenhofF 
& Heidt ti2002i hereafter MHOl) found bulge shape indices of 
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0.8 < n < 6.7 for their sample of 40 bright galaxies, with only 
one object having a bulge with « < 1 and two objects with 
n > 4. 

The fits were judged to be reliable {x^, < 2) in 1 17 (of 174) 
galaxies. The average percentage errors for the fitted parame- 
ters of bulge and disk are less than 20%, consistent with MHOl 
and MACH03. We have further restricted the final database 
to contain only those galaxies with apparent optical diameter 
a>V (from NED'). Thus, the final PP sample consists of 108 
galaxies, 2 of which do not have well-defined morphological 
types, and to which we have assigned T=5 (following LEDA, 
the Extragalactic Database based in Lyon). The best fit param- 
eters resulting from our 2D fits for these 108 disk galaxies are 
reported in Table 1 . Distances to individual galaxies are deter- 
mined from the Hubble law (with Ho - 70 km s"' Mpc"') after 
correcting the heliocentric radial velocities listed in Paper I to 
the CMB reference frame according to the prescription of the 
Third Reference Catalogue (de Vaucouleurs et al. 1991 here- 
after RC3). The full table is only available in electronic form, 
please see the footnote to the Abstract, here we reproduce it 
partly to illustrate its contents. Table 1 also contains additional 
information from RC3 as follows: 

Column 1: UGC galaxy number, if not differently speci- 
fied; 

Column 2: Distance (Mpc); 
Column 3: Hubble type T; 
Column 4: bar class ; 
Column 5: global of fit; 

Column 6: bulge shape index n (« = for pure disks); 
Column 7: apparent bulge ellipticity e^; 
Columns 8-9: bulge effective radius re and uncertainty (arc- 
sec); 

Columns 10-11: bulge effective surface brightness, reduced to 

face-on, He and uncertainty (mag arcsec"^); 

Column 12: total bulge apparent magnitude m^ (mag); 

Column 13: disk inclination / (degrees); 

Columns 14-15: disk scalelength and uncertainty (arcsec); 

Columns 16-17: central face-on disk surface brightness fiQ (and 

uncertainty mag arcsec"^); 

Column 18: total disk apparent magnitude m^ (mag); 
Column 19: bulge-to-disk luminosity ratio B/D; 
Column 20: total galaxy apparent magnitude m (mag). 

2.2. Photometric corrections 

The magnitudes have been corrected for Galactic extinction ac- 
cording to the estimates given by Burstein & Heiles ( 1984 1 as 
reported by NED, and assuming the extinction curve of Cardelli 
et al. ([1989 1. Using the recipe of Frogel et al. ( I1978t . we have 
applied corrections for redshift dimming. 

While internal dust extinction may affect the estimates of 
structural parameters even in the NIR (e.g., Peletier & Willner 
[1992 1, the corrections are small and subject to some uncer- 
tainty. An empirical analysis of the photometric corrections for 

' The NASA/IPAC Extragalactic Database (NED) is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space Administration. 



dust extinction in this sample has been caiTied out by Moriondo 
et al. ( 1998b hereafter Paper II). To do this, they rely on corre- 
lations with inclinations of the observed disk and bulge struc- 
tural parameters and colors, and find that the structural param- 
eters of the bulge are not affected by dust extinction, while 
those of the disk may be - especially the disk scalelength - but 
with large scatter However, the analysis Paper II is based on 
the assumption that dust is homogeneously distributed within 
a thin disk coincident with the mid-plane of the stellar distri- 
bution, which in general is not true (e.g. Bianchi et al. 119961 
2000 FeiTara et aI. fT999) . Moreover, the behavior of the kine- 
matic scalelength as a function of the axial ratio is complex 
(Giovanelli & Havnes 120021 see also Paper II). Therefore, for 
the subsequent analysis, the structural parameters (r^, fi^, 
r^, m^, B/D and m) are corrected to face-on assuming full 
transparency, which implies no corrections to scalelengths, to- 
tal magnitudes, and bulge-to-disk luminosity ratio B/D. 

All the above corrections have already been applied to the 
data reported in Table 1 . 

2.3. Comparison witii previous work 

Our sample contains five galaxies (UGC 463, UGC 927, 
UGC 1089, UGC 12378, UGC 12527) in common with 
MACH03 (H band). In all cases, the two sets of disk parameters 
agree quite well: r^(ours)/r£/(theirs) = 0.95 + 0.05 and //o(ours)- 
yL(o(theirs) = 0.16 + 0.06 H-mag arcsec"^. For UGC 927, they 
classify their profile as a "truncated disk", but it appears very 
different from our major-axis cut (all our fits can be obtained in 
graphic form by request from the authors). They find n = 2.7, 
the largest in their sample, while we find n - A; the bulge sur- 
face brightnesses differ by about ~ 1 mag. This may be an ex- 
ample of the pitfalls of ID fitting since this galaxy has twisted 
isophotes, and marked spiral structure even in the NIR. Such 
morphological features tend to be suppressed by our fitting al- 
gorithm, and the greater number of degrees of freedom in 2D 
image fitting is clearly an advantage. 

In general, we find larger n values than those by MACH03: 
the mean difference An = 1.1 + 0.75. We also find larger 
bulge r^, but this is expected because of the different n's. This 
may be another advantage of our fitting algorithm, as it takes 
into account the non-spherical shape of the bulge (e.g., MGH). 
Bulges are rarely spherical (Kent 1988 1, and not considering 
the possible bulge ellipticity would result in a steeper bulge 
outer profile, and thus a lower n. 

3. Results 

A priori the parent sample suffers from a known dearth of early- 
type spkals (e.g., r < 2, see Paper I). Indeed our present sam- 
ple contains only 6 galaxies with type Sab or earlier, with a 
sample median of r = 5 (Sc). Therefore, we have added to the 
subsequent analysis the sample of 14 early-type spirals (Sa's) 
studied by MGH; the same fitting algorithm was used to de- 
compose the images.^ The combined sample comprises 122 

- A correction of 0.25 mag has been made to convert the K-band 
quantities in that paper to the H band used here (e.g., Giovanardi & 
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galaxies, now with 21 of type Sab or earlier, and with median 
T = 4(Sbc). 

Throughout the analysis, we will make use of mean effec- 
tive surface brightnesses (SFBs), calculated as the average flux 
within a circular area of radius r^. The disk scalelengths are 
converted to effective ones by multiplying by a factor of 
1.6784 appropriate for n - \ (see MGH), and we define the 
mean effective SFB < /if > and < yuf > for bulge and disk re- 
spectively. In the following, to avoid confusion, the bulge scale- 
length will be indicated as rf . Mean SFBs < fie > were chosen, 
rather than the central /iq (Khosroshahi et al. 2000a 2000b|, or 
IJe - the SFB at (e.g., Graham 2001 MACH03) - because 
of the strong functional dependence on n of the latter param- 
eters. If a given bulge is fit with a generalized exponential, jj^ 
becomes fainter as « of the fit becomes larger, because of the 
mathematical form of the generalized exponential. Higher n fits 
are more strongly peaked and consequently have a higher yUo, 
but are more extended, hence a larger r^, a fainter yu^, and a 
higher bulge luminosity (cf.. Fig. 2 in MGH). In practice, 
when the same bulge is fit with n varying between 1 and 4, the 
resulting /j.^ changes typically by more than 3 magarcsec"^, 
while < fig > also becomes fainter but hy ^ I magarcsec"^. 
The adoption of < yu^. > also conforms to the definitions of RC2 
and RC3, and to all the lore of the Fundamental Plane (FP). 

3.1. Pure disk galaxies 

Only two of the galaxies (UGC 975 and UGC 1579) are best fit 
with pure disks {B/D = 0); bulges are found in all the others, 
although in some cases as faint as 16.9 H mag, with B/D ra- 
tios of 1%. UGC 975 is classified as a generic spiral ("S"), as it 
shows on blue images a tight flocculent pattern and high surface 
brightness. UGC 1579, an almost face-on SB(s)d, shows an ex- 
tended spiral structure and a large central region of relatively 
high surface brightness. They are among the faintest objects in 
the sample, but still galaxies with M ~ -22 H mag, and are 
not low surface brightness galaxies. Pure disks are rare in our 
sample, although it comprises 50 galaxies with T>5. The selec- 
tion criteria adopted tend to avoid very late types of low surface 
brightness; for this reason, our two pure disks are unlikely to be 
representative of their class (cf. Gavazzi et al. .2000J . 

3.2. Possible biases 

In principle, the performance of our decomposition procedure 
could depend on the size and quality of the image. We have 
performed a simple check for such biases by examining the dis- 
tributions of apparent image parameters versus the bulge shape 
index n, as shown in Fig. |2l optical properties are taken from 
RC3. The only evidence is that: (/) the two pure-disk (n = 0) 
galaxies are faint, and (//) the Sa's from MGH are substantially 
brighter than the PP set. For the rest, no obvious significant 
trends are detectable. In particular, we observe no trend for 
faint, small, and highly inclined galaxies, that is those observed 
with lower relative resolution, to be better fit with any particular 

Hunt [1996 (. The distance scale has also been transformed to the Ho 
(70 kms"' Mpc"') adopted for this study. 
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Fig. 2. Distributions of H-band apparent magnitude m (top 
panel), optical axial ratio b/a (middle), and apparent optical 
size a (lower) as a function of bulge shape index n. Galaxies 
best fitted with n = 1,2,3,4 are shown by x, open squares, 
open triangles, and filled circles, respectively. The early-type 
spirals from MGH are also included and shown as open circles 
independently of their n. Data for pure-disk galaxies are shown 
by -H at n = 0. Large open squares correspond to mean values, 
and the error bars to Icr standard deviation. 

« values. This is important for the reliability of any assertions 
regarding trends with bulge shape index. 

3.3. Hubble type and bulge shape n 

The upper panel of Fig. |3] shows the distribution of Hubble 
types T for our combined sample. The PP and MGH samples 
have different shadings and the paucity of early-type spirals in 
the original PP sample is evident. 

The lower panel shows the means of the bulge shape index 
n averaged over different Hubble types. Here and in the fol- 
lowing we have grouped Hubble types into four bins: T < 2, 
2 < r < 4, 4 < r < 6, and r > 7. It can be seen from the 
figure that n - 1 bulges tend to be more common in later-type 
spirals, and that early types (T < 2) tend to have bulges with 
larger n, significantly different from n = I. This result agrees 
with other authors (Andredakis et al. ll9951 Graham 2001 > who 
found that n tends to increase with decreasing Hubble type. 
Such a trend could also explain the inconsistency with the 
results of MACH03 who found no evidence for bulges with 
n > 2, but whose sample contains only two galaxies earlier 
than Sb (r < 3). 

In what follows, we will examine trends of structural pa- 
rameters both with the bulge index n and with Hubble type 
T . The justification for this is that n and T are only weakly 
correlated (see Fig|3l and previous work). We want to inves- 
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Fig. 3. Distribution of Hubble type T (upper panel) and bulge 
shape index n versus T (lower panel). In both panels, as in all 
subsequent figures, the galaxies from MGH are included. In 
the histogram, the shading for the PP sample is right slanted 
and for the MGH sample is left slanted. In the lower panel, 
horizontal error bars indicate the range of Hubble type over 
which the mean n is calculated. The vertical bars are Icr stan- 
dard deviation. The horizontal dotted line is placed at n = 1 for 
convenience. The number of objects in each T bin is printed 
over its data point. 
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per panel), mean effective surface brightnesses (middle), and 
absolute magnitude (lower). The right panels show the same 
quantities for the disks. Symbols as in Fig.|2] 



tigate trends with bulge shape as well as with Hubble type 
since bulges with different n may be intrinsically different (e.g., 
Carollo lT999l Carollo et al. l200il . independently of the Hubble 
type of their host galaxies. Therefore, with the aim of singling 
out the more significant parameter, bulge shape n or Hubble 
type T, we will analyze trends of photometric parameters with 
both. 

3.4. Bulge/disk parameters versus bulge shape n 

The behavior of the structural parameters of bulges and disks 
as a function of the shape index n of the bulge is illustrated in 
Fig El Inspection of the figure shows that neither bulge nor disk 
< > vary significantly with n. Obviously, when the bulge /i^ 
or yUo are used (not shown), there is a large variation which is a 
result of the functional form of the fitting law. The increase of 
central bulge surface brightness /io with n is ~ 6 mag arcsec"^ 
in agreement with Khosroshahi et al. (2000a I. 

rf , instead, increases significantly with n: the left panel of 
Fig- S shows that n - \ bulges have smaller scalelengths than 
bulges with higher tr', while there is no dependence of disk 

^ This effect is significant at > 99.9% level (U test). However, the 
significance is only obtained for n = \ bulges; the differences among 
rf at other n values are not significant. 



scalelengths on the n of the bulge. The increase of with n is 
consistent with earlier work on ellipticals (Caon et al. I1992> . 
on early-type spirals (Khosroshahi et al. ,2000a i. and on the 
later-type spirals in the de Jong sample (de Jong & van der 
Kruit 1994; Graham 2001 ). Such a trend was shown to be real, 
rather than a result of the mutual dependence among n, rf , and 
He of the bulge fitting law, by Trujillo et al. ( 1200 U . Moreover, 
since the properties of n = 1 bulges are shown to differ with 
high statistical significance from their higher n counterparts, it 
vindicates our method which relies on integer values of n. 

The lower panels of Fig. 0] shows and as a func- 
tion of n. In agreement with earlier work (Andredakis et al. 
[T9951 Graham 2001., MHOl; Trujillo et al. 2002), bulges with 
larger n tend to be more luminous, although with some scatter 
Again, we find that n - \ bulges are significantly different from 
bulges with higher n, being less luminous with a significance 
> 99.9%. On the other hand, appears to be independent of 
n, and significantly less variable in general: the standard devi- 
ations of range from 1.2 to 1.6 mag, while those of 
are considerably smaller, 0.5 to 0.8 mag. The SFBs do not vary 
with n, and cannot be the cause of the dependence of on 
n. Exponential bulges are smaller and less luminous than their 
higher n counterparts, but roughly the same in terms of mean 
surface brightness. None of the disk properties (< /if >, rf , 
M^) depends on the n value of the bulge. 
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Fig. 5. Bulge and disk structural parameters versus Hubble type 
T. The left panels show bulge effective scale lengths (upper 
panel), mean effective surface brightnesses (middle), and abso- 
lute magnitude (lower). The right panels show the same quan- 
tities for the disks. Symbols as in Fig.|2] 

3.5. Bulge/disk parameters versus Hubble type 

Here we examine the same parameters of the previous section, 
but now consider their variation with Hubble type T. The rele- 
vant graphs are presented in Fig.|5] 

The bulge < //f > significantly brighter with decreasing 
Hubble type; the effect is formally highly significant^, a new 
result made possible by the inclusion of the MGH Sa's, and the 
consequently more ample range in T, as well as by our use of 
mean SFBs < ji^ >, rather than /y^ . Although de Jong J1996> 
found a similar trend, all bulges were fit with simple exponen- 
tials; it was therefore not clear if the effect was real for early- 
type bulges which are better fit with higher n. No such trend 
was noted by MHOl, although as they point out, their sample 
extends only to T = 5. Also the disk < fif > gets brighter with 
decreasing Hubble type, although the trend is much weaker 
and not readily discernible in Fig. |5f . This also is a new re- 
sult, which again arises from the large range of Hubble types 
in our combined sample. 

Regarding both bulge and disk effective scalelengths, 
the variations with T are small. Although very late-type disks 
appear to be smaller than their early-type counterparts, both 

Early spiral types {T < 2), intermediate (2 < T < 4) and late 
(4 < r < 6), and very late (T > 7) types differ at a significance 
level of > 99.9%, with the exception of late and very late types whose 
difference is significant at ~ 94%. 

' The increase in disk SFB is formally significant only for early 
(T < 2) and intermediate (2 < T < 4) types, relative to very late types 
(T > 7); the first is > 99.9% significant, while the latter is significant 
only at a ~ 95% level. 



Fig. 6. Bulge versus disk mean surface brightnesses < fie >, 
separately for galaxies with different best n values of the bulge. 
The correlation coefficients are shown below the values of n. 
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Fig. 7. Bulge versus disk effective scalelengths rf, rf, sepa- 
rately for galaxies with different best n values of the bulge. The 
correlation coefficients are shown below the values of «. 



bulge and disk scalelengths are relatively constant with Hubble 
type, unlike the variation of bulge rf with n discussed in the 
previous section. For a given Hubble type, the range spanned 
by disk properties (M^, < fif >, ) is smaller than that of the 
bulge (M*, < juf >, rf ). 
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Fig. 8. Bulge versus disk absolute magnitudes M , M , sepa- 
rately for galaxies with different best n values of the bulge. The 
correlation coefficients are shown below the values of n. 

The lower panels of Fig.|5]plot the absolute bulge and disk 
magnitudes and as a function of T. There is a clear 
trend for the bulge luminosity to decrease going from early to 
very late-type spirals^. A similar result was also obtained by 
MHOl . Absolute disk magnitudes do not vary significantly with 
T, except that very late types T>7 are less luminous than all 
earlier type s^. 

3.6. Bulge parameters versus disk parameters 

Figures |6lEl and |8l illustrate the correlations, separated by n, 
between the parameters of a bulge and those of its host disk. For 
all parameters, there is a clear trend for n - I bulges and disks 
to be tightly correlated. For n - 2, the correlation is somewhat 
degraded, and by the time n - 3,4, the correlations have all but 
disappeared. The tendency is evident for mean SFBs (Fig.|6j, 
scalelengths (Fig.Q, and luminosities (Fig.|8}. 

While this is a new result, one of the main ones of this 
paper, previous work hinted at something similar. With ID 
bulge-hdisk decompositions of R-band images, and assigning 
n - 1 or 4, Andredakis & Sanders ( 1994) found a correlation 
between bulge and disk absolute magnitudes, with the correla- 
tions being tighter for systems with n = 1 than for n - 4 bulges. 
MHOl found weak correlations between //f and //q, and with rf 
and disk scalelength , but did not separate or bin the n val- 
ues; thus the trends evident in our data were diluted in their 
continuous distribution. 

* The significance of the trend is >99.9% for all bins, except that 
T<2 and slightly later types are not formally distinguishable. 

' Significance is >99% except for T<2 which has ~95% signifi- 
cance). 



Fig. 9. Bulge mean effective surface brightness < yuf > versus 
effective radius rf in kpc (upper panel). Disk mean effective 
surface brightnesses < yuf > versus effective radius rf (lower 
panel). Symbols as in Fig. |2 The best-fit regression (OLS bi- 
sector) for all data is shown in both panels. 



Correlations between bulge and disk radii (Fig.0 are not 
new, and have been used to argue for secular evolution as the 
primary mechanism of bulge formation (e.g., Courteau et al. 
I1996> . Independently of the ratio of radii, which will be dis- 
cussed in the following section, our data show that the corre- 
lation between radii is best for n = 1, but degrades steadily 
toward higher n values. By the time n = 4, the correlation has 
disappeared completely. If a sample does not contain higher n 
bulges because of a dearth of early-type spirals, the bulge/disk 
radii will appear correlated (e.g., MACH03). Our result also 
agrees with Graham (2001 ) who found a constant ratio of bulge 
and disk radii only up to n ~ 2. Khosroshahi et al. (EOOdat and 
Scodeggio et al. <2002> also report the existence of such a cor- 
relation although they found large scatter; again, this is likely 
due to the dilution of the correlations when all n values are 
combined together as done by these authors. 

3.7. Bulge/disk scaling relations 

Figure|5]shows mean surface brightness < fig > versus effective 
radius, separately for bulges and disks. We find both bulge and 
disk < fie > to be well correlated with their effective radii, in 
the sense that larger bulges and disks tend to be fainter in SFB. 
This is well known; such correlations have long been found for 
spiral bulges (Kent '1985 , Kodaira et al. '1986', Andredakis et 
al. .1995j . and constitute the "Kormendy" relation for ellipti- 
cals (Kormendv l 19771 Hoessel & Schneider ll985l Djorgovski 
& Davis "1989't. More recent studies (MHOl; Khosroshahi et 
al. i2000aj also find a similar correlation but with larger scatter. 



L. K. Hunt et al.: On the origin of bulges 



9 



cq 






12 3 4 



5 
T 



Fig. 10. Bulge and disk effective scalelength ratios rf /rf (up- 
per panels) and bulge-to-disk luminosity ratios (lower) versus 
n (left panels) and Hubble type T (right). Symbols as in Fig.|2 



Bulges and disks appear to share the same FP, but occupy dis- 
tinct, though overlapping, regions of it (Burstein et al. .1997j . 
Combining bulges and disks together, we obtain the regres- 
sion (shown in the figure) < >- 16.81 + 3.21 logr^, 
that is re oc 7^**^**, which is consistent with canonical esti- 
mates of the FP relation (e.g., Kormendy & Djorgovski 1989 1: 
re °^ /-0-83±o.08 fjjggg ^j-g j^Qth Substantially steeper than the 
ti-end reported in MHOl (r^ oc '^). 

SFBs plotted against scalelengths form an almost face-on 
view of the FP (e.g., Kormendy & Djorgovski I1989> . It can 
be seen in Fig. |9] that for bulges (top panel) the main scatter 
in the bulge relation is introduced by n = 1 values (x, r - 
0.30), while the best correlation is for n - 4 (filled circles, 
r - 0.80). As in Carollo ( 1999 1, for a given radius, n - I bulges 
tend to be fainter than their higher n counterparts. As a whole, 
disks are significantly better correlated (r = 0.71) than bulges 
(r = 0.53 when all n are included). Although this projection of 
the FP is not ideal for such considerations, these results imply 
that residuals relative to the FP may be correlated with shape 
parameters (e.g., Hjorth & Madsen [19951 Prugniel & Simien 
[T9971 Khosroshahi et al. 2000b i. 

Since scalelength ratios constitute the observational cor- 
nerstone for secular evolution scenarios (Courteau et al. 1 19961 
MACH03), we shall discuss this parameter in detail for our 
data set. Figure [TUl shows B/D and rf/rf as a function of n 
(left panels) and T (right panels). While there is a clear trend 
for rf /rf to increase with n, there is only a very weak trend 
for smaller rf/rf with increasing T. For n - \ bulges, we 
find < rf/rf >= 0.085 + 0.037 (51 galaxies); for n = 2 
< rf/rf >= 0.15 +0.10 (31); for n = 3 < rf /rf 0.21 +0.16 



(18); for n = 4 < rf /rf 0.32 + 0.32 {I9f. The size ratio in- 
creases from n = 1 to « = 4 by a factor of ~ 4. The mean value 
of /rf we find for n = 2, 3 is consistent with the MACH03 
value of 0.22 + 0.09 for Type-I profiles, but this value is more 
than twice as large as for the n - \ bulges in our sample. 
Independently of Type I vs. Type II profiles, we surmise that 
much of the scatter in the ,rf relation is due to the mixing of 
diff'erent n bulges. As was seen in Section IT6l the correlation 
between bulge and disk scalelengths is best for n - I bulges, 
then gradually decays to a non-correlation for n - A (see Fig. 
0. In addition to the trend of increasing < rf /rf > with n, the 
dispersion also increases with n: n — I bulge/disk size ratios 
have a scatter ~ 9 times smaller than that for n = 4. 

While < rf/rf > depends on n, it is rather independent 
of Hubble type. The means < /rf > for different T bins 
are statistically indistinguishable, although very late-type spi- 
rals (r > 7) do tend to have < r^ jr^ > smaller than very early 
types {T <2): 0.12+0.07 (10 galaxies) vs. 0.19+0.13 (19 galax- 
ies). MACH03 also find a slight dependence of < /rf > on 
morphology (~ 0.20 - 0.24 from late- to early-type spirals), 
although the trend is very weak, and not confirmed by Graham 
(2001), nor with formal significance by our data. Hence, al- 
though variations of /rf with T are diflicult to detect, this 
ratio changes significantly with the shape of the bulge, imply- 
ing that the spiral sequence cannot be scale free (c.f., Courteau 
et al.[l996 MACH03). This point will be discussed further in 
the context of secular evolution scenarios in Section^] 

As expected, we find that B/D increases with n and 
decreases with T, in agreement with previous work (e.g., 
Andredakis et al. 119951 Khosroshahi et al. I2000al Graham 



120011 MHOl). B/D is a criterion for the Hubble classification, 
so the correlation with T is obvious. The dependence of B/D 
on n is similar, although of lower significance, than its depen- 
dence on r^. In any case, the spread of B/D for a given T bin is 
large, approaching two orders of magnitude, and the looseness 
of B/D vs. T is rather surprising. It is clear that the Hubble 
classification reflects complex, multivariate processes, which 
are not easily quantified. 

B/D is larger for early-type spirals because of a brighter 
< fie >' since bulge rf , and the disk parameters do not change 
significantly with T. This is a confirmation of the "iceberg" 
scenario proposed by Graham (I2001t . On the other hand, B/D 
is larger for larger n bulges because rf is larger, since the disk 
parameters and < yuf > are roughly constant. When galaxies 
are divided according to the shape of their bulges rather than 
their Hubble type, the spiral sequence is seen to be not scale 
free. It appears therefore that the iceberg scenario proposed by 
Graham (Ej 



cannot be only a function of surface bright- 



* The difference in effective scalelength ratio among n values is 
>99.9% significant for n = 1 and larger ns, and >99% for n = 2 
and n = 4. n = 2 and n = 3 ratios do not differ significantly, nor do 
n = 3 and n = 4. 

' For our data, the significance of the trend with T is formally 99% 
or greater between all type bins, except for early type spirals (T < 2) 
and the next bin (2< T <4), where it is slightly lower (~97%). n = 
1 bulges have a significantly (>99.9%) lower B/D ratio than higher 
n values, but the differences among other n values are not formally 
significant. 
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Fig. 11. Upper panels: galaxy absolute magnitudes as a func- 
tion of bulge shape n (left panel) and of Hubble type T (right); 
symbols as in Fig.|2 Lower panels: mean difference in absolute 
magnitude between strongly barred (SB) and unbarred (SA) 
galaxies vs. bulge shape index n (left panel) and of Hubble type 
T (right). 



ness. The shape of the bulge is a main factor for its increasing 
prominence with decreasing T; this is true not only because of 
the variation of , but also because of the variation in bulge 
central surface brightness /i^ due to the functional form of the 
Sersic law. 



3.8. Luminosities 

Bearing in mind that our galaxies are rather bright objects, with 
total galaxy luminosity M"" (H band) i - 22, roughly equiva- 
lent to M'"' (B band) i - 18, the sample spans roughly 4 mag- 
nitudes, providing some dynamic range for statistical analysis. 
The upper panels of Fig. shows the total absolute galaxy 
magnitude (sum of bulge and disk components) vs. bulge shape 
n (left panel) and Hubble type T (right). As known for some 
time (e.g., Roberts & Havnes rT994l l. there is a significant trend 
of luminosity with Hubble type, in that early-type spirals tend 
to be more luminous than late types'". We find only a weak 
trend of luminosity with «, since galaxies with n - I bulges 
tend to be less luminous than higher n galaxies at the ~ 95% 
level, and that the two pure disks are among the least luminous 
objects in the sample. In any case, the range in luminosities for 
a given Hubble type bin is large (~ 3 mag), implying that other 



In our sample, this trend is >99% significant for all bins relative 
T >1\ the difference in luminosity between 2 < T < 4 and 4 < T < 6 
is significant at > 95%, although T < 2 and 2 < T < 6 do not have 
significantly diff^erent luminosities. 



Fig. 12. Bulge (solid circles) and disk (open circles) 
absolute magnitudes as a function of the total absolute magni- 
tude of the galaxy. Hence, there are two data points for each 
galaxy. 



variables are at work in addition to T, e.g., DDO classification 
(van den Bergh [T960t . 

That earlier Hubble types tend to be more luminous is ex- 
plained by the analogous trend in bulge luminosity only for 
T < 1, since for these types disks are of roughly the same lu- 
minosity. However, for T > 1 disks tend to be less luminous 
than those of earlier types, so that the lower total luminosities 
are due to a lower luminosity of both components, not just of 
the bulge. 

and are plotted against M"" in Fig.[^ The average 
luminosities of bulge and disk increase with the luminosity of 
the parent galaxy, independently of the shape n of the bulge. 
The slope of the correlation is unity for disks and a rather steep 
2.5 for bulges. For disks associated with n - I bulges, the cor- 
relation is extremely tight since ~ M'°' for these systems. 
In addition, ~ M'"' for the less luminous disk galaxies, 
while ~ for the most luminous ones. In other words, in 
the most luminous galaxies, the bulge-to-disk luminosity ratio 
B/D~ 1. 

3.8.1 . Influence of bars 

We have attempted to assess the influence, if any, of bars 
on galaxy luminosity, and trends with Hubble type and bulge 
shape index. In the lower panels of Fig. ^2 we plot the dif- 
ference between the mean absolute luminosities < Ms b > and 
< Ms A > of the sample galaxies classified in RC3 as SB (strong 
bar) and S A (no bar), respectively; the left panel shows this dif- 
ferences plotted vs. n, and the right one vs. T. Inspection of the 
figure clearly shows that strongly barred galaxies with expo- 
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nential bulges are significantly more luminous than their un- 
barred counterparts. Despite the trend of increasing luminosity 
with n shown in Fig. ^2 the 6 strongly barred n — \ galaxies 
are the most luminous subset in our sample: < Ms sin = 1) > = 
-24.5 + 0.9. There is also a continuous trend toward smaller 
differences < Ms b > - < Ms a > with increasing n. However, 
the same is not true for trends with Hubble type (right panel); 
we find no evidence that a strong bar makes galaxies of a given 
type more (or less) luminous. 



4. Discussion 

First, we summarize the main results obtained in the previous 
section. 



- Bulge effective scalelengths rf increase with increasing n; 
n - 1 bulges are smaller than those with n > 1 (Fig. 
Neither bulge and disk < ju^ > nor disk rf vary with n. 
Bulge luminosity M^ increases with n, while disk luminos- 
ity is independent of n (Fig.|3. In particular, exponen- 
tial bulges are smaller and less luminous than those with 
n > I. 

- Bulge < juf > gets brighter as T decreases (Fig. Disk 
< fj.^ > shows a similar, although weaker, trend. Bulge 
luminosity increases with decreasing T (Fig. |5jl, but 
only disks with T>7 are significantly less luminous than 
their early-type counterparts. Neither bulge nor disk scale 
lengths rf , vary significantly with Hubble type (Fig.|5|i- 

- When bulge parameters (< fj.^ >, , M*) are compared 
with disk ones (< /if >, rf , M^), in all cases they are 
tightly correlated for n - \ bulges (Figs.|6lE]|H}. The cor- 
relations gradually worsen with increasing n such that n - A 
bulges appear virtually independent of their disks. 

- The scatter in the correlation between bulge < yuf > and 
rf (Fig. |9} is shown to depend on bulge shape «; the two 
parameters are tightly correlated in n = 4 bulges (r = 0.8), 
and increasingly less so as n decreases. Disk < yuf > and 
rf are well correlated (r = 0.7). 

- Luminosity ratios B/Daie smaller for n = 1 bulges than for 
n > 1, and increase with decreasing T, but with a surpris- 
ingly large spread (Fig. llO> . The value of B/D is governed 
by the bulge luminosity, since disks have the same average 
luminosity independently of their shape index n. Bulge-to- 
disk size ratios rf/rf are also smaller for n - \ bulges, 
but are independent of Hubble type. The value of rf / rf for 
n = \ bulges is 4 times smaller than that for n - 4, with a 
spread which is 9 times smaller 

- As already known, early-type spirals tend to be more lumi- 
nous than late types (Fig. II U: except for T >7, this is due 
to the increasing contribution of the bulge. 

- Strongly barred galaxies with exponential bulges are more 
luminous than their unbarred counterparts. The difference 
tends to vanish with increasing n. There is no analogous 
trend with T. 



4.1. Bulge formation and evolution 

The better correlations of bulge and disk properties for expo- 
nential bulges suggest that n - \ bulges and their disks are 
closely related. The tightness of the relation worsens as n in- 
creases, so that when n = 4, bulges appear independent of their 
disks. In contrast, the (Kormendy) relation between bulge SFB 
and scalelength is very tight for n - A bulges, and increasingly 
noisy for lower n values. This implies that « = 4 bulges are in 
some sense more structurally defined than those with lower n. 

In the Introduction, we sketched the competing scenarios 
for bulge formation; here we use the results of our study in 
an attempt to distinguish among them. In exclusively "early" 
bulge formation scenarios, in which the bulge formed before 
the disk, it is difficult to understand the close relation we ob- 
serve between exponential bulges and their disks. This rela- 
tion and its gradual degradation as n increases is not easily 
explained by either the monolithic collapse scenario or CDM 
merging of primordial disks. If initial conditions were every- 
thing, such as in a monolithic collapse, bulge and disk prop- 
erties should be statistically unrelated. The same would be 
true if stable disk growth depended on satellite accretion af- 
ter bulge formation from merging disks, such as in the CDM 
merging scenarios. It may be that such processes operate for 
some (higher n) bulges, but they are not a good explanation for 
the global properties of low n ones. Nor can they explain the 
gradual worsening of bulge/disk correlations with n. 

A close relation between bulges and disks is, however, fore- 
seen by the inside-out formation scheme of van den Bosch 
( I1997> . In this scenario, the collapsing halo forms unstable 
disks which, through bar instabilities, rapidly become "bulges". 
A disk is able to form only after the bulge becomes large 
enough to stabilize the ensuing bar. Nevertheless, in this sce- 
nario, it is difficult to explain the existence of luminous galax- 
ies with very low BID ratios. If the presence of a disk de- 
pends on the prior buildup of a sufficiently luminous (massive) 
bulge, then we would not expect to find galaxies with luminous 
disks and relatively small bulges. Indeed, the existence of pure 
disk galaxies (~2% of our sample, but see Gavazzi et al. l2000> 
would appear to contradict this scenario. 

A more natural explanation of our observations is the sec- 
ular evolution scheme in which n - \ bulges are relatively 
"young" structures, having formed "recently" from their disks 
through internal gravitational instabilities. This process must 
be dissipative because of phase-space arguments (e.g., Wyse 
.1998) . but may not exclusively involve bar-driven gas inflow 
(e.g., Pf ennige r & Norman [T9^ C ombes et al. [19901 Friedli 
& Benz II 99*31 Norman et al.'1996). The evolutionary scheme 
proposed by Zhang in a series of papers (Zhang 1996 1998. 
1 19991 1 predicts that "as a galaxy evolves through the spiral- 
induced collective dissipation process, its surface density distri- 
bution (i.e., disk plus bulge) will become more and more cen- 
trally concentrated, together with the buildup of an extended 
outer envelope". Moreover, in a given galaxy, the observed 
spread in the angular momentum distribution reflects the degree 
to which secular evolutionary processes have been at work. 
Since stable stellar systems require a greater degree of veloc- 
ity anisotropy as n decreases (Ciotti &. Lanzoni .l997J . n = 1 
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bulges must be kinematically less isotropic than higher n ones. 
If this anisotropic component is related to tangential motion, 
and thus to angular momentum, such bulges must be -in the 
Zhang scenario- younger. 

Let us therefore hypothesize that exponential bulges are 
indeed younger than those of higher «, independently of the 
Hubble type of the galaxy in which they reside. In other words, 
the first spheroids to be formed from their disks would be 
bulges with n - I. Such a hypothesis is supported by the blue 
V-H colors of « = 1 bulges, relative to those with n - 4 (CaroUo 
et al. 2001 1. It would also be consistent with the good correla- 
tion of exponential bulge properties with those of their disks as 
we have seen in Sect. 13.61 If indeed the bulge n increases with 
time, as predicted by the Zhang models, then our observation 
that rf and increase with n implies that they also increase 
with time. The good n = 4 correlation of < fi^ > with rf would 
then be a consequence of their older age and better defined "re- 
laxed" (virialized?) structure, and their consequent better con- 
sistency with the FP. The degradation of the correlations be- 
tween bulge and disk parameters with increasing n would also 
imply that larger-n bulges are increasingly more distant in time 
from their disks. In such a scheme, the bulge and disk branches 
in Fig.^lcould be interpreted as the "evolutionary tracks" of 
these components. 

4.2. The Hubble sequence 

Appealing as this picture may be, the trends with Hubble type 
still need to be explained. Consistently with previous work, 
we find that higher n values are more common in early-type 
spirals. Conversely, all (our sample) galaxies with T >7 have 
« = 1 bulges. Therefore, there is some indication that n and 
Hubble type are (albeit weakly) correlated, and may both be 
associated with age or evolutionary state. Since the degree of 
evolution is expected to be partly determined by galaxy mass 
(Zhang [1999 1, early type spirals which are also more lumi- 
nous (massive) should in some sense be more evolved than 
late types. That surface brightness increases with decreasing 
Hubble type for bulges and to a lesser extent for disks is con- 
sistent with this expectation. For a given (constant) rate of star 
formation from the formation epoch to the present day, and a 
sufficient gas reservoir, we would expect higher stellar surface 
densities in galaxies which have been forming stars for longer 
If, on the other hand, star-formation rates were to depend on 
(gas-Hstellar) mass, then we would expect higher mass systems 
to have higher stellar densities. Either way, it is difficult to rec- 
oncile the higher stellar surface densities in early type spirals 
(Fig. 13 with an earlier evolutionary state (or younger age); 
they are probably more evolved (i.e., older) systems. Early-type 
spirals also tend to have lower gas fractions (e.g., Roberts & 
Haynes[l994i than later types, presumably a result of having 
converted the available gas into stars. 

If late-type spirals (n = 1 bulge) are younger than early 
types (higher n), the question of bars and their effect on evo- 
lution must be addressed, since in local samples of galax- 
ies, >75% of very late-type spirals {T >7) are barred (Ho et 
al. [19971 Hunt & Malkan J999j . The excess of bars in very 



late-type systems may contribute to thek morphological evo- 
lution as proposed by many groups (e.g.. Combes et al. 119901 
Pfenniger & FriedH 1991 FriedH & Benz 1993 Hasan et al. 
119931 Norman et al. 1996 1. Numerical models show that bulge- 
like structures with an exponential n - \ profile can be pro- 
duced from bar-like perturbations that induce inward radial in- 
flow of disk matter (Pfenniger & Friedli 119911 Zhang & Wyse 
2000 1. Our finding that strongly barred n - \ bulges tend to be 
more luminous than unbarred ones could reflect such evolution. 

Indeed, this bar-driven evolutionary scheme has been ad- 
vocated by Courteau et al. (I1996> and MACH03 to explain the 
relative constancy of B ID size ratios they find in (the Type I 
profiles of) their late-type spiral sample. However, they claim 
that the ratio rf/rf is constant for all spiral types, with the 
implication that the Hubble sequence is scale free. One of the 
main results of our study is that several of the trends observed 
are clearly detected as a function of bulge shape «, rather than 
T\ we interpret this as due to the complex interplay of the spi- 
ral classification criteria. With a sample containing few early- 
type {T < 2) spirals, and thus few high-n bulges, several of 
the trends we found would not be detectable. Hence, we find a 
wider range of n, in agreement with other groups (Andredakis 
et al. 1995] MHOl; Graham ^001 1, and also a wider range of 
rf/rf . We also find that the bulge-to-disk size ratio changes 
significantly with n. While such variations with n indicate that 
the Hubble sequence is not scale free, they do not necessarily 
disprove the secular evolution scenario. However, the scenario 
must be expanded to include the genesis of bulges with differ- 
ent (higher) n's, together with their distinct properties. 

Bar-driven inflow as the unique mechanism of secular evo- 
lution may not be sufficient to explain the full range of prop- 
erties in the Hubble sequence of spirals. The analytical mod- 
els of Zhang & Wyse (JOOOl based on viscous evolution of 
bars find no galaxies with BID ~ 1; their bulges are not 
sufficiently luminous compared to their disks. Moreover, the 
secular-evolution simulations of Bouwens et al. ( 1999 1 which 
rely on similar assumptions also produce bulges which are too 
small and too blue. The collective stellar dissipation scheme 
proposed by Zhang (1999') may also be insufficient. While 
qualitatively consistent with our findings, this scheme produces 
eff'ects (e.g., buildup of central concentration) which are much 
smaller than those we observe, perhaps because the eff'ects of 
gas are not included. 

Evidently, some other mechanism is needed, such as merg- 
ers (e.g., Baugh et al. 1996) or accretion episodes (e.g., 
Steinmetz & Navarro .2002 1. However, even mergers in hier- 
archical clustering scenarios are not able to generate the range 
of bulge/disk properties in our sample. Scannapieco & Tissera 
(2003 1 use hydrodynamical simulations to assess the effects of 
mergers on disk galaxy structure. While they are able to repro- 
duce the bulge-to-disk size relation for the n < 2 bulges in our 
sample and in that of MACH03, the tight correlation of rf /rf 
for « = 1 bulges is not explained. Moreover, they fail to predict 
the strong increase of rf with n (Fig.|3]and other work), since 
their models predict the opposite trend. Nevertheless, they suc- 
cessfully reproduce the BID ratio as a function of bulge shape, 
although their n values tend to cluster around 1, unlike the 
early-type spirals in our sample which tend to be best fit by 
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larger n values (Fig.|3- On the other hand Aguerri et al. (12001 1 
are able to reproduce n - A bulges with satellite accretion, but 
only with satellites as massive as the bulge itself. 

4.3. Conclusions 

The hypothesis that emerges from our study is essentially the 
following: 

- Disks are formed first, during an initial collapse, and have 
properties which depend on the mass and angular momen- 
tum of the dark-matter (DM) halo. 

- Bulges are then born first as n = 1 systems, through secular 
transfer of disk material to the central regions, probably by 
bar instabilities. 

- With an evolutionary rate that depends on initial conditions 
(gas-i-DM mass, angular momentum), n of the bulge in- 
creases, with consequent growth of rf , M^, and B/D, per- 
haps by means of the collective dissipation mechanism de- 
scribed by Zhang ( 1999 1, together with the effects of gas 
transfer 

- As n gets larger and the bulge builds up, Hubble type gets 
earlier, with consequent brightening (SFB) of bulges and 
(less so) disks, and an increase of M*, B/D. 

- Since n - \ bulges are relatively unevolved systems, in 
a state of formation at the present epoch, they are struc- 
turally not as relaxed as larger n bulges, and present a large 
scatter in the < yuf > - log rf (Kormendy) relation. Their 
stellar luminosity is enhanced by the presence of a strong 
bar driving the evolution. The n - \ bulges that we see 
today are those systems which have had long evolutionary 
times, probably because of their lower luminosity (mass). 
The evolutionary state of the disk would also be a func- 
tion of n, since the entire system would be expected to have 
evolved faster when n is larger 

- Mergers of bulge/disk systems and accretion of gas proba- 
bly affect the buildup of the bulge, but are probably impor- 
tant only for the creation of high-n bulges. 

- The correlation between "evolutionary state indicators" n 
and T is not expected to be perfect because of the complex 
nature of spiral classification, and because of variations in 
galaxy mass, angular momentum, bar strength, gas fraction, 
dark halo merging history, and perhaps other parameters. 

- The gradual growth of a bulge as n increases is implied by 
the continuity we observe in bulge properties as a function 
of n. All trends are tight for n - 1, and steadily worsen 
with increasing n, except for the Kormendy relation which 
is best followed by structurally evolved systems (n = 4). 

One prediction of such a scenario is a certain spread in 
stellar age for high-n bulges. Detailed simulations are neces- 
sary to assess whether this prediction is in conflict with the 
level of synchronicity in the star-formation history observed in 
spheroids (Bower et al. 1992). We also would expect that ex- 
ponential bulges, prevalent in late-type spirals and closely con- 
nected with their disks, would be more dominated by rotation 
than their higher n counterparts. Low-luminosity spheroids are 
known to have a higher degree of rotational support than more 



luminous ones (Davies et al. 1983), so it may be that this is true 
also for n = 1 bulges relative to those of higher n. 

Clearly more modelling is needed to understand the full 
range in properties along the Hubble sequence, and what they 
imply for bulge formation. However, secular evolution appears 
to be a strong driver of the Hubble sequence. Our observations 
are compatible with the hypothesis that n and T are both rough 
measures of bulge age, and that bulges are, at least initially 
(« - 1), formed secularly from disks. 
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20.72 


0.64 


14.86 


11 


4.6 


0.2 


17.72 


0.13 


12.43 


0.11 


12.32 


833 


74 


6 


SB 


1.2 


3 


0.02 


1.1 


0.8 


18.19 


0.87 


14.83 


46 


12.6 


0.6 


18.17 


0.06 


10.68 


0.02 


10.66 


835 


98 


6 


none 


1.3 


1 


0.29 


0.5 


0.2 


17.91 


0.83 


16.84 


46 


10.7 


0.4 


18.71 


0.04 


11.56 


0.01 


11.55 


886 


84 


5 


SAB 


1.5 


1 


0.33 


2.3 


0.1 


17.75 


0.04 


13.22 


49 


13.9 


0.3 


18.30 


0.03 


10.58 


0.09 


10.49 


909 


72 


7 


none 


1.1 


1 


0.26 


2.3 


0.3 


19.32 


0.15 


14.78 


47 


11.0 


0.4 


18.77 


0.05 


11.57 


0.05 


11.52 


919 


78 


2 


none 


1.2 


1 


0.37 


2.2 


0.1 


17.08 


0.07 


12.69 


69 


16.1 


0.2 


17.64 


0.02 


9.61 


0.06 


9.55 


927 


86 


4 


none 


1.1 


4 


0.23 


2.6 


0.4 


19.41 


0.19 


13.96 


55 


12.2 


0.2 


18.42 


0.03 


11.00 


0.07 


10.93 


940 


100 


5 


SA 


1.0 


1 


0.15 


1.6 


0.3 


19.30 


0.19 


15.52 


49 


10.4 


0.4 


19.05 


0.05 


11.97 


0.04 


11.93 


962 


72 


4 


SB 


1.6 


1 


0.10 


3.3 


0.1 


16.09 


0.01 


10.82 


27 


26.0 


0.3 


17.87 


0.01 


8.80 


0.16 


8.64 


975 


70 


5 


none 


1.4 





0.00 


0.0 


0.0 


0.00 


0.00 


0.00 


52 


6.3 


0.1 


17.75 


0.02 


11.75 


0.00 


11.75 


1005 


59 


6 


none 


1.4 


4 


0.21 


32.9 


10.2 


22.54 


0.36 


11.57 


51 


18.6 


0.9 


18.79 


0.12 


10.44 


0.35 


10.12 


1013 


74 


3 


SB 


1.5 


3 


0.43 


13.6 


0.2 


18.60 


0.02 


9.70 


79 


78.1 


15.2 


22.17 


0.18 


10.71 


2.53 


9.34 


1033 


58 


6 


none 


1.6 


2 


0.78 


25.4 


1.7 


21.11 


0.12 


11.05 


84 


5.8 


0.4 


18.67 


0.13 


12.87 


5.33 


10.87 


1034 


74 


4 


SB 


1.2 


1 


0.28 


1.2 


0.1 


17.03 


0.13 


13.97 


64 


12.8 


0.2 


17.79 


0.02 


10.27 


0.03 


10.23 


1089 


71 


6 


none 


1.0 


4 


0.15 


9.2 


0.7 


20.40 


0.07 


12.19 


33 


10.2 


0.1 


17.89 


0.04 


10.85 


0.29 


10.58 


1100 


65 


3 


SAB 


1.8 


3 


0.34 


5.7 


0.2 


18.12 


0.04 


11.12 


75 


29.5 


1.3 


19.70 


0.06 


10.36 


0.50 


9.92 


1131 


72 


8 


SAB 


0.9 


1 


0.05 


2.1 


0.6 


19.97 


0.28 


15.68 


19 


9.9 


1.1 


18.83 


0.13 


11.86 


0.03 


11.83 


1210 


77 


3 


none 


1.7 


4 


0.19 


2.0 


1.0 


18.68 


0.70 


13.78 


37 


6.4 


0.5 


17.24 


0.17 


11.22 


0.09 


11.12 


1238 


64 


3 


none 


1.2 


2 


0.22 


6.3 


0.1 


17.86 


0.02 


10.83 


47 


24.0 


2.8 


20.09 


0.16 


11.20 


1.41 


10.25 


1250 


53 


-2 


SA 


1.1 


3 


0.53 


12.0 


0.3 


18.64 


0.02 


10.00 


74 


11.5 


0.3 


18.61 


0.07 


11.31 


3.35 


9.72 


1259 


74 


5 


SAB 


1.8 


1 


0.23 


2.1 


0.1 


17.02 


0.04 


12.72 


46 


14.3 


0.3 


17.72 


0.03 


9.95 


0.08 


9.86 


1291 


49 


5 


SAB 


1.2 


1 


0.38 


3.6 


0.7 


19.83 


0.16 


14.33 


53 


15.7 


0.8 


18.78 


0.07 


10.81 


0.04 


10.77 


1349 


87 


6 


none 


1.2 


2 


0.01 


8.4 


2.4 


19.68 


0.11 


12.03 


45 


14.8 


4.8 


19.38 


0.62 


11.53 


0.63 


11.00 


1376 


80 


4 


SB 


1.1 


2 


0.41 


3.0 


0.4 


19.21 


0.15 


13.76 


56 


17.0 


0.7 


18.78 


0.05 


10.63 


0.06 


10.57 


1395 


74 


3 


SA 


1.5 


4 


0.22 


1.3 


0.1 


16.87 


0.06 


12.87 


48 


12.7 


0.2 


18.23 


0.02 


10.73 


0.14 


10.58 


1451 


70 


3 


SB 


1.5 


1 


0.01 


1.9 


0.1 


15.46 


0.01 


11.37 


66 


11.6 


0.2 


17.72 


0.02 


10.40 


0.41 


10.03 


1564 


90 


4 


SB 


1.3 


3 


0.51 


2.9 


0.3 


18.46 


0.12 


12.89 


63 


7.1 


0.2 


17.50 


0.06 


11.25 


0.22 


11.03 


1579 


49 


7 


SB 


1.1 





0.00 


0.0 


0.0 


0.00 


0.00 


0.00 


31 


10.7 


0.2 


18.79 


0.01 


11.64 


0.00 


11.64 
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i 


rd 




fJ-o 






B/D 


m 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


1582 


69 


7 


none 


0.9 


1 


0.00 


3.1 


1.0 


19.41 


0.19 


14.28 


76 


14.3 


2.2 


19.46 


0.16 


11.69 


0.09 


11.60 


1633 


60 


5 


SAB 


1.5 


1 


0.42 


2.0 


0.1 


16.82 


0.07 


12.66 


62 


12.1 


0.1 


16.95 


0.02 


9.54 


0.06 


9.48 


1688 


73 


2 


SB 


1.8 


3 


0.24 


3.0 


0.5 


17.46 


0.21 


11.80 


42 


12.1 


0.3 


16.94 


0.04 


9.52 


0.12 


9.40 


1695 


64 


3 


none 


1.4 


1 


0.11 


3.4 


0.2 


16.80 


0.04 


11.45 


82 


12.4 


0.4 


17.95 


0.06 


10.48 


0.41 


10.11 


1768 


56 


4 


SAB 


1.1 


3 


0.35 


6.5 


0.6 


18.49 


0.09 


11.19 


51 


16.5 


0.3 


17.21 


0.03 


9.13 


0.15 


8.98 


1792 


71 


5 


SAB 


1.1 


4 


0.00 


5.0 


2.0 


19.55 


0.46 


12.65 


54 


15.9 


0.8 


18.55 


0.09 


10.55 


0.14 


10.40 


1835 


72 


5 


none 


2.0 


2 


0.03 


3.6 


0.5 


17.95 


0.12 


12.13 


52 


16.4 


1.1 


18.22 


0.09 


10.15 


0.16 


9.99 


1887 


79 


5 


SA 


1.0 


1 


0.34 


0.9 


0.4 


16.23 


0.65 


13.66 


50 


12.0 


0.3 


17.75 


0.03 


10.35 


0.05 


10.30 


1937 


59 


6 


none 


1.3 


4 


0.17 


14.6 


3.2 


19.95 


0.17 


10.74 


57 


15.7 


0.5 


17.75 


0.09 


9.78 


0.41 


9.40 


2042 


61 


5 


SA 


1.6 


4 


0.19 


9.3 


0.5 


18.15 


0.05 


9.92 


37 


6.7 


0.1 


16.38 


0.04 


10.25 


1.36 


9.32 


2094 


73 


5 


SB 


1.0 


4 


0.02 


5.7 


1.4 


20.49 


0.28 


13.33 


11 


15.5 


0.9 


19.04 


0.09 


11.10 


0.13 


10.97 


2109 


52 


7 


SAB 


1.1 


1 


0.10 


1.9 


0.2 


19.04 


0.13 


14.94 


27 


18.1 


0.5 


18.88 


0.03 


10.59 


0.02 


10.57 


2156 


64 


5 


SA 


1.1 


1 


0.21 


1.5 


0.2 


17.72 


0.16 


14.11 


39 


11.4 


0.3 


18.01 


0.03 


10.72 


0.04 


10.68 


2178 


56 


1 


SB 


1.6 


4 


0.11 


2.9 


0.5 


17.12 


0.20 


11.42 


49 


11.1 


0.3 


17.53 


0.07 


10.30 


0.36 


9.97 


2185 


62 


6 


none 


1.1 


1 


0.15 


0.9 


0.1 


17.30 


0.17 


14.73 


75 


16.9 


0.3 


18.39 


0.02 


10.26 


0.02 


10.25 


2204 


65 


5 


SAB 


1.1 


1 


0.04 


2.9 


0.2 


18.21 


0.09 


13.24 


16 


13.0 


1.0 


18.46 


0.10 


10.89 


0.11 


10.77 


2223 


71 


6 


none 


1.4 


2 


0.03 


1.8 


0.7 


20.15 


0.71 


15.79 


74 


11.4 


0.9 


18.85 


0.10 


11.57 


0.02 


11.55 


2241 


97 


4 


SA 


0.9 


1 


0.20 


4.8 


0.3 


17.97 


0.04 


11.90 


51 


10.7 


1.7 


18.44 


0.28 


11.30 


0.58 


10.80 


2264 


91 


5 


SB 


1.1 


1 


0.20 


0.7 


0.1 


17.54 


0.09 


15.50 


54 


11.5 


0.2 


18.48 


0.02 


11.19 


0.02 


11.17 


2353 


51 


3 


SAB 


1.3 


1 


0.09 


1.5 


0.1 


15.45 


0.04 


11.81 


45 


12.0 


0.2 


17.07 


0.02 


9.68 


0.14 


9.53 


2368 


51 


3 


SB 


1.9 


2 


0.17 


2.4 


0.3 


16.79 


0.10 


11.86 


76 


23.1 


1.5 


19.04 


0.07 


10.22 


0.22 


10.01 


2435 


69 


6 


SA 


1.1 


3 


0.29 


18.6 


4.1 


21.69 


0.18 


12.10 


46 


36.2 


3.4 


20.08 


0.15 


10.30 


0.19 


10.11 


2548 


68 


3 


SA 


1.4 


3 


0.06 


3.2 


0.2 


17.01 


0.07 


11.26 


37 


16.6 


0.5 


17.95 


0.05 


9.85 


0.27 


9.59 


2586 


63 


3 


none 


1.8 


4 


0.50 


24.0 


2.2 


20.85 


0.12 


10.56 


62 


4.8 


0.1 


15.96 


0.04 


10.54 


0.99 


9.80 


2604 


65 


5 


SAB 


1.0 


1 


0.30 


2.3 


0.3 


18.91 


0.10 


14.41 


47 


9.2 


0.3 


18.02 


0.05 


11.20 


0.05 


11.15 


2617 


69 


7 


SAB 


1.1 


1 


0.38 


1.9 


0.3 


18.17 


0.15 


14.04 


61 


12.4 


0.5 


18.51 


0.06 


11.05 


0.06 


10.98 


2627 


60 


5 


SA 


1.7 


1 


0.00 


2.8 


0.1 


17.55 


0.04 


12.59 





14.2 


0.3 


17.91 


0.06 


10.16 


0.11 


10.05 


2655 


88 


7 


SAB 


1.2 


2 


0.38 


12.7 


4.3 


20.16 


0.15 


11.60 


53 


17.4 


6.1 


19.60 


0.62 


11.41 


0.84 


10.75 
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Bar 




n 




re 




l^e 
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ra 










B/D 


m 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


2658 


44 


3 


SAB 


1.5 


1 


0.25 


1.1 


0.1 


17.17 


0.17 


14.22 


43 


8.4 


0.3 


17.77 


0.04 


11.15 


0.06 


11.09 


2810 


93 


5 


SA 


1.4 


1 


0.04 


1.1 


0.2 


17.42 


0.18 


14.51 


41 


10.8 


0.4 


18.22 


0.04 


11.06 


0.04 


11.01 


2885 


83 


5 


SA 


1.5 


3 


0.01 


4.3 


0.2 


17.36 


0.06 


10.95 


68 


34.7 


0.9 


18.72 


0.03 


9.03 


0.17 


8.85 


2931 


82 


5 


SA 


1.1 


2 


0.26 


3.6 


1.1 


19.59 


0.22 


13.79 


43 


8.3 


0.7 


18.52 


0.20 


11.92 


0.18 


11.74 


11931 


66 


4 


SB 


1.2 


1 


0.09 


1.2 


0.2 


16.52 


0.27 


13.43 


25 


10.6 


0.5 


17.43 


0.05 


10.32 


0.06 


10.26 


12103 


102 


3 


SB 


1.3 


1 


0.21 


1.4 


0.1 


15.43 


0.09 


12.08 


45 


7.6 


0.2 


16.00 


0.04 


9.61 


0.10 


9.50 


12108 


80 


3 


none 


1.4 


2 


0.14 


2.3 


0.3 


18.49 


0.12 


13.63 


31 


7.2 


0.2 


17.11 


0.05 


10.84 


0.08 


10.76 


12137 


67 


4 


SAB 


1.1 


1 


0.08 


3.6 


0.2 


18.53 


0.07 


13.08 


23 


26.0 


1.4 


19.04 


0.06 


9.96 


0.06 


9.90 


12173 


68 


5 


SAB 


1.6 


1 


0.14 


1.4 


0.3 


17.49 


0.25 


14.09 


49 


17.3 


0.5 


17.99 


0.03 


9.81 


0.02 


9.79 


12230 


92 


4 


none 


1.3 


1 


0.00 


1.3 


0.2 


17.32 


0.13 


14.00 


58 


8.8 


0.2 


17.28 


0.04 


10.57 


0.04 


10.52 


12286 


104 


5 


SA 


1.4 


3 


0.00 


1.5 


0.3 


17.19 


0.30 


13.13 





8.3 


0.3 


17.46 


0.09 


10.87 


0.13 


10.74 


12378 


89 


7 


none 


1.1 


1 


0.19 


1.2 


0.2 


17.35 


0.13 


14.32 


55 


9.8 


0.2 


17.75 


0.04 


10.80 


0.04 


10.75 


12379 


89 


4 


none 


1.1 


4 


0.09 


8.3 


0.3 


18.39 


0.04 


10.41 


25 


18.6 


0.5 


18.32 


0.05 


9.98 


0.67 


9.42 


12486 


71 


4 


SB 


1.1 


2 


0.05 


4.2 


0.1 


16.88 


0.02 


10.75 


64 


22.7 


1.2 


19.19 


0.07 


10.41 


0.73 


9.81 


12527 


61 


6 


none 


1.6 


1 


0.34 


1.8 


0.2 


17.09 


0.13 


13.16 


50 


11.6 


0.3 


17.52 


0.03 


10.21 


0.07 


10.14 


12539 


54 


3 


SA 


1.2 


1 


0.23 


1.6 


0.1 


16.95 


0.03 


13.25 


65 


11.1 


0.1 


17.47 


0.02 


10.24 


0.06 


10.17 


12586 


61 


3 


none 


1.0 


4 


0.09 


9.8 


0.9 


18.59 


0.09 


10.24 


33 


10.7 


0.6 


18.30 


0.19 


11.15 


2.32 


9.85 


12598 


50 


5 


none 


1.5 


1 


0.15 


1.2 


0.1 


16.39 


0.09 


13.26 


59 


7.9 


0.1 


16.51 


0.02 


10.02 


0.05 


9.97 


12618 


73 


-2 


SB 


1.0 


2 


0.00 


3.5 


0.2 


16.36 


0.05 


10.64 


36 


11.4 


1.5 


18.50 


0.25 


11.22 


1.70 


10.14 


12666 


71 


6 


none 


1.2 


1 


0.33 


5.5 


0.8 


19.70 


0.08 


13.28 


74 


28.7 


6.6 


20.17 


0.18 


10.88 


0.11 


10.77 


12667 


54 


6 


none 


1.4 


1 


0.38 


1.3 


0.3 


18.23 


0.28 


14.92 


53 


14.0 


0.4 


18.16 


0.03 


10.44 


0.02 


10.42 


12780 


74 


4 


SAB 


1.7 


1 


0.06 


2.0 


0.1 


16.20 


0.03 


11.98 


40 


16.0 


0.2 


17.35 


0.02 


9.33 


0.09 


9.24 


12821 


96 


6 


none 


1.1 


3 


0.08 


4.4 


1.2 


19.59 


0.28 


13.14 


24 


9.7 


0.5 


18.06 


0.11 


11.13 


0.16 


10.97 


12915 


62 


5 


none 


1.6 


1 


0.32 


0.6 


0.3 


16.04 


0.79 


14.43 


65 


6.4 


0.1 


16.55 


0.04 


10.51 


0.03 


10.48 



